How to Transfer between Arbitrary n-Qubit Quantum States by Coherent Control 
and Simplest Switchable Noise on a Single Qubit 
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We explore reachable sets of open n-qubit quantum systems, the coherent parts of which are 
under full unitary control and that have just one qubit whose Markovian noise amplitude can be 
modulated in time such as to provide an additional degree of incoherent control. In particular, 
adding bang-bang control of amplitude damping noise (non-unital) allows the dynamic system to 
act transitively on the entire set of density operators. This means one can transform any initial 
quantum state into any desired target state. Adding switchable bit-flip noise (unital), on the other 
hand, suffices to explore all states majorised by the initial state. We have extended our open-loop 
optimal control algorithm (dynamo package) by such degrees of incoherent control so that these 
unprecedented reachable sets can systematically be exploited in experiments. As illustrated for 
an ion trap experimental setting, open-loop control with noise switching can accomplish all state 
transfers one can get by the more complicated measurement-based closed-loop feedback schemes. 

PACS numbers: 03.67.-a, 03.67.Lx, 03.65.Yz, 03.67.Pp; 89.70.+C 
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Recently, dissipation has been exploited for quantum 
state engineering p], H[ so that evolution under constant 
noise leads to long-lived entangled fixed-point states. 
Lloyd and Viola Q showed that closed-loop feedback 
from one resettable ancilla qubit suffices to simulate any 
quantum dynamics of open systems. Both concepts were 
used to combine coherent dynamics with optical pumping 
on an ancilla qubit for dissipative preparation of entan- 
gled states [ij or quantum maps Q. — Besides, full con- 
trol over the Kraus operators @ or the environment 0| 
allows for interconverting arbitrary quantum states. 

Manipulating quantum systems with high precision is 
paramount to exploring their properties for pioneering 
experiments, and also in view of new technologies Q. 
Therefore it is most desirable to extend the current tool- 
box of optimal control Q by systematically incorporating 
dissipative control parameters. 

Here we prove that it suffices to include as a new con- 
trol parameter a single bang-bang switchable Markovian 
noise amplitude for one qubit (no ancilla) into an oth- 
erwise noiseless and coherently controllable network to 
increase the power of the dynamic system such that any 
target state can be reached from any initial state. Thus 
we extend our optimal-control platform dynamo 
by controls over Markovian noise sources. To illustrate 
possible applications, we demonstrate the initialisation 
step [Tl| of quantum computing, i.e. the transfer from the 
thermal state to the pure state |00 . . . 0) (as well as the 
opposite process) , the interconversion of random pairs of 
mixed states, and finally the noise-driven generation of 
maximally entangled states. 

Theory. We consider the quantum Markovian master 
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equation of an n-qubit system as a bilinear control sys- 
tem (£): 

p(t) = -{iH u + T)p(t) and p(0) = Po (1) 



with H u := Ho + V , '/,! /)//, 
evolution Hamiltonian Hq, the 



comprising the free- 
control Hamiltonians 



Hj switched by piecewise constant control amplitudes 
Uj(t) € R and H u as the corresponding commutator su- 
peroperator. Take T to be of Lindblad form 



r(p) 



E {VtpVi - \{V}Vtp + P V}V,)) , (2) 



where now 7e(i) <E [0,7*] with 7* > will be used as 
additional piecewise constant control parameters. 

In the sequel we will consider mostly systems with a 
single Lindblad generator. In the non-unital case it is the 
Lindblad generator for amplitude damping, V a , while in 
the unital case it is the one for bit flip, V b , defined as 



V a 



If"" 1 ® |0)(1| 



and V b :=l\ 



.<j x /2. (3) 



Here we follow the Lie-algebraic setting along the lines 
of [lMll- As in we say the control system on n 
qubits meets the condition for (weak) Hamiltonian con- 
trollability if the Lie closure under commutation is 

(iH , iHj I j = 1, ... , m) Ue = su(N) with N := 2™. (4) 

Now the reachable set teacf) s (/?o) is defined as the set of 
all states p(T) with T > that can be reached from po 
following the dynamics of (£). If Eqn. §4$ holds, without 
relaxation one can steer from any initial state po to any 
other state /Qtarget with the same eigenvalues. In other 
words, for 7 = the control system (E) acts transitively 
on the unitary orbit U(p Q ) := {U p U^ \ U € SU(N)} of 
the respective initial state po. This holds for any po in 
the set of all density operators, termed pos 1 henceforth. 
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Under coherent control and constant noise (7 > non- 
switchable) it is difficult to give precise reachable sets 
for general n-qubit systems that satisfy Eqn. (j4]) only 
upon including the system Hamiltonian (Hq). Based on 
seminal work by Uhlmann majorisation criteria 

that are powerful if Hq is not needed to meet Eqn. (|4]) 
[20I . I2H now just give upper bounds to reachable sets by 
inclusions. But with increasing number of qubits n, these 
inclusions become increasingly inaccurate and have to be 
replaced by Lie-semigroup methods as described in [22| . 

In the presence of bang-bang switchable relaxation on 
a single qubit in an n-qubit system, here we show that 
the situation improves significantly and one obtains the 
following results, both proven in the Supplement [23j|: 

Theorem 1. Let T, a be an n-qubit bilinear control sys- 
tem as in Eqn. ([T]) satisfying Eqn. ([¥]) for 7 = 0. Suppose 
the n th qubit (say) undergoes (non-unital) amplitude- 
damping relaxation, the noise amplitude of which can be 
switched in time between two values as j(t) G {0, 7*} with 
7* > 0. If the free evolution Hamiltonian Hq is diagonal 
(e.g., Ising-ZZ type), and if there are no further sources 
of decoherence, then the system S a acts transitively on 
the set of all density operators pos 1 : 



teacf) Ea (/9 ) = pos 1 for all p € pos 1 , 

where the closure is understood as the limit 7*T 



00. 



Theorem 2. Let be an n-qubit bilinear control system 
as in Eqn. (JTJ) satisfying Eqn. Q (7 = 0) now with the 
n th qubit (say) undergoing (unital) bit- flip relaxation with 
switchable noise amplitude j(t) £ {0,7*}. If the free evo- 
lution Hamiltonian Hq is diagonal (e.g., Ising-ZZ type), 
and if there are no further sources of decoherence, then 
in the limit 7»T — > 00 the reachable set to explores all 
density operators majorised by the initial state po, i.e. 



ceocf) Eb (p ) = {p G pos 1 I p -< pa} for any p Q G pos 1 . (6) 

The conditions for the drift Hamiltonian Hq in the the- 
orems above can be relaxed. The details are given in the 
Supplement |23|, E], along with the proofs. 

The scenarios of Eqn. (|3]) can be generalised to the 

Lindblad generator Vg := with 9 G [0,1]. If 

9 7^ 1/2 the noise qubit has a unique fixed point Pao(9). 
Comparing this to the canonical density operator of a 
qubit with energy level splitting A, relaxation by Vg gives 
the same fixed point as equilibrating the system via the 
noisy qubit with a local heat bath of inverse temperature 
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artanh (8(6)) with 8(9) 



e 2 



using the shorthand 9 := 1 — 6. As limiting cases, pure 
amplitude damping is brought about by a bath with zero 
temperature (i.e. 9 — 0), while pure bit-flip corresponds 
to the high-temperature limit T — > 00 (i.e. 9 -|). 

While a single-qubit system with unitary control and 
bang-bang switchable noise generator Vg can clearly be 



steered (asymptotically) to any state with purity less or 
equal to the larger of the purities of po and Poo{9), the 
situation for n > 2 qubits is more involved: using coher- 
ent control, relaxation of a diagonal initial state can be 
limited to a single pair of eigenvalues at a time if all the 
remaining ones can be arranged in pairs (pa , pjj ) , each 
satisfying 



2 /e 2 < pu/pjj < 9 2 /9 2 



(7) 



Yet Eqn. ([7]) poses no restriction in the important task 
of cooling: starting from the maximally mixed state, op- 
timal control protocols with period-wise relaxation by Vg 
interspersed with unitary permutation of diagonal den- 
sity operator elements is more general than the partner- 
pairing approach [24j to algorithmic cooling with bias 
8(9) and < 9 < 1/2. This type of algorithmic cooling 
proceeds also just on the diagonal elements of the density 
operator, but it involves no transfers limited to a single 
pair of eigenvalues (details in the Supplement (23l [B]). 

Exploring Model Systems. To challenge our extended 
optimal control algorithm, we first consider two exam- 
ples of state transfer where the target states are on the 
boundary of the respective reachable sets of the initial 
states, in other words, they can only be reached asymp- 
totically (7*T — > 00). To illustrate Theorems Q] and [2] we 
then demonstrate noise-driven transfer (i) between ran- 
dom pairs of states under controlled amplitude damping 
noise and (ii) between random pairs of states satisfying 
Ptarget -< Po under controlled bit-flip noise in Examples 
3 and 4. We finish with entanglement generation in an 
experimental trapped ion system in Example 5. 

In Examples 1-4, our system is an n-qubit chain with 
uniform Ising-ZZ nearest-neighbour couplings given by 

Ho '■= ttJ ^TJfc o^z^ 1 ^ and piecewise constant x and y 
controls (that need not be bounded) on each qubit locally, 
so the control systems satisfy Eqn. (j4]). We add control- 
lable noise (either amplitude-damping or bit flip) with 
amplitude j(t) G [0,7*] acting on one terminal qubit. In 
all the examples we set 7* = 5 J. 

Example 1. Here, as for initialising a quantum com- 
puter the task is to turn the thermal initial state 
Pth '■= -prl into the pure target state |00 . . . 0) by unitary 
control and controlled amplitude damping. For n qubits, 
the task can be accomplished in an n-step protocol: let 
the noise act on each qubit q for the time T q to popu- 
late the state |0)(0|, and permute the qubits between the 
steps. A linear chain requires 57J" =1 (g — 1) = (™) nearest- 
neighbour swaps. Since all the intermediate states are 
diagonal, the swaps can be replaced with i-swaps, each 
taking a time of -j under the Ising-ZZ coupling. The 
residual Frobenius error Sf is minimised when all the T q 
are equal, giving 



(e) = l-2(l- 



(1- 



¥ 2 ) 



(8) 



where e := e t*T n /n an( j ._ YJ^r g . Linearizing this 
expression and adding the time for the i-swaps, the total 
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Figure 1. Transfer from the thermal state p t h = gl to the 
zero-state P|ooo) m a 3-qubit Ising-ZZ chain with controlled 
amplitude-damping noise on qubit one as in Example 1. 
(a) Quality versus total sequence duration T. The dashed line 
gives the upper bound from Eqn. @, and the dots (red circles 
for averages) individual optimization runs with random initial 
sequences. Noise amplitudes were initialised in three distinct 
blocks of equal duration to help the optimisation towards an 
economic solution, (b) Evolution of the eigenvalues under 
the best of the T = 6/ J solutions. This sequence (c) shows 
three relaxative periods with maximal noise amplitude 7 a i for 
transforming eigenvalues, while unitary actions governed by 
(u xv ,u yv ) mainly take place in the intervals between. Each 
purely unitary segment is of the approximate duration 1/J, 
corresponding to the duration of a single i-swap. 

duration T a of this simple protocol as a function of Sp 
amounts (in first order) to 

r -(;)*+* h (^F)- w 

Fig. [T] demonstrates that optimal control can outper- 
form this simple scheme by parallelising part of the uni- 
tary transfer with the amplitude-damping driven 'cool- 
ing' steps. Interestingly, the initialisation task can still 
be accomplished to a good approximation when unavoid- 
able constant dephasing noise on all the three qubits is 
added, as shown in the Supplement (23. ID]. 

Example 2. In turn, consider 'erasing' the pure initial 
state 1 00 . . . 0) to the thermal state p t h- Under controlled 
amplitude damping, this can be accomplished exactly, 
each round splitting the populations in half with a to- 
tal time of T' a = Q)^ + ^M 2 )- However, with bit-flip 
noise this transfer can only be obtained asymptotically. 
One may use a similar n-step protocol as in the previous 
example, this time approximately erasing each qubit to 
a state proportional to 1. Again, optimal control greatly 
outperforms this simple scheme. Results and details are 
shown in the Supplement [2^, [D] . 
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(a) (b) (c) 




wall time [s] time [1/J] time [1/J] 



Figure 2. (a) Quality vs. computation time for state trans- 
fer between pairs (po, ptarget) of random 3-qubit states using 
controlled amplitude-damping noise (solid) in addition to lo- 
cal unitary control. Same for random pairs (po,Ptar g et) with 
Ptargct -< po under controlled bit flip noise (dashed). In both 
cases (a) shows the median of 9 optimisation runs for each of 
the eight random state pairs. Representative examples of evo- 
lution of the eigenvalues for an amplitude-damping transfer 
(b) and for a bit-flip transfer (c). In the former case, a typical 
feature is the initial zeroing of the smaller half of the eigenval- 
ues while the larger half are re-distributed among themselves. 
Only at the very end are the smaller eigenvalues resurrected. 

Example 3. We illustrate transitivity under controlled 
amplitude damping on one qubit plus general unitary 
control by transfers between pairs of random 3-qubit den- 
sity operators. Fig.^a) shows the algorithm to converge 
well to Sf — 10 -4 . As shown in Fig. [H^b), the best se- 
quences seem to zero the smaller half of the eigenvalues as 
soon as possible just to revive them in the very end after 
the larger half has been balanced among themselves. 

Example 4. Similarly allowing for controlled bit-flip 
noise on one qubit plus general unitary control, we ad- 
dress the transfer between arbitrary pairs of 3-qubit den- 
sity operators with ptarget ^ Po- Fig. IH a ) again illus- 
trates the good convergence of the algorithm. Many of 
these profiles exhibit "terraces" which could indicate local 
quasi-optima. The unital case may be harder to optimise 
in general: (1) the majorisation condition entails that a 
suboptimal transfer made early in the sequence cannot 
be outbalanced later in the control sequence; it can only 
be mended in a next iteration; (2) the necessity for si- 
multaneous decoupling (like Trotterisation in the proof 
of Thm. [5]) adds to the hardness of the optimisation. 

Example 5. The final example addresses entanglement 
generation in a system similar to the one in It con- 
sists of four trapped ion qubits coherently controlled by 
lasers. On top of individual local z-controls (u z \, ■ • ■ , u z ^) 
on each qubit, one can pulse on all the qubits simultane- 
ously by the joint x and y-controls F v := | X^=i a vj with 
v = x, y as well as by the quadratic terms F 2 := (F y ) 2 . 
All the control amplitudes are expressed as multiples of 
an interaction strength a. In contrast to where the 
protocol resorts to an ancilla qubit to be added (fol- 
lowing Q) for a measurement-based circuit on the 4+1 
system, here we do without the ancilla qubit by making 
just the terminal qubit subject to controlled amplitude- 
damping noise with strength 7 al , to drive the system 
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Figure 3. State transfer from the thermal state to the four- 
qubit GHZ state in the ion-trap system of Example 5 similar 
to 0| . By controlled 'pumping' (amplitude damping) on one 
qubit, one can do without closed-loop measurement-based cir- 
cuits involving an additional ancilla qubit as required in [3,3- 
Our sequence (a) drives the system to the state (c), which dif- 
fers from the target state IGHZ4) by an error of 5f — 5 ■ 10~ 3 . 
The time evolution of the eigenvalues (b) illustrates parallel 
action on all the eigenvalues under the sequence. 



from the thermal initial state p t h '■— to the pure en- 
tangled target state |GHZ 4 ) = 75GOOOO) + 1 1111)). As 
shown in Fig. [3J the optimised controls use the noise with 
maximal amplitude over its entire duration interrupted 
just by two short periods of purely unitary control. 

Discussion. By unitary controllability, we may diago- 
nalise the initial and the target states. So transferring a 
diagonal initial state into a diagonal target state can be 
considered as the normal form of the state-transfer prob- 
lem. It can be treated analytically, because it is easy to 
separate dissipation-driven changes of eigenvalues from 
unitary coherent actions of permuting eigenvalues and 
decoupling drift Hamiltonians. Now the difference be- 
tween optimising amplitude-damping non-unital transfer 
(as in Thm. [IJ and bit-flip unital transfer (as in Thm. ^ 
becomes evident: In the non-unital case, transitive ac- 
tion on the set of all density operators clearly helps to 
escape from suboptimal intermediate control sequences 
during the optimisation. Yet in the unital case, the ma- 
jorisation condition ptargot -< p(t) -< po for all < t < T 
and the boundary conditions p(0) = po, p(T) = ^target 
(at worst for 7*7 — > 00) explain potential algorithmic 
traps: one may easily arrive at an intermediate state 
pm{t) -< po that comes closer to the target state, but 
will never reach it as it fails to meet the reachability con- 
dition ptargot -< Pm(t); see the Supplement [23T.IE] for how 
to avoid this. 

One may contrast our method with the closed-loop 
control method in Q originally designed for quantum- 
map synthesis using projective measurement of a coupled 



resettable ancilla qubit plus full unitary control to enact 
arbitrary quantum operations (including state transfers) , 
with Markovian evolution as the infinitesimal limit. Ap- 
plied to state transfer, the present method instead relies 
on a switchable local Markovian noise source and requires 
no measurement nor an ancilla [25| . 

Conclusions and Outlook. We have proven that by 
adding as a new control parameter bang-bang switch- 
able Markovian noise on just one system qubit, an other- 
wise coherently controllable n-qubit network can explore 
unprecedented reachable sets: in the case of amplitude- 
damping noise (or any noise process in its unitary equiv- 
alence class, with compatible drift) one can convert any 
initial state po into any target state ptarget, while un- 
der switchable bit-flip noise (or any noise process uni- 
tarily equivalent) one can transfer any po into any tar- 
get ptargct -< Po majorised by the initial state. These 
results have been further generalised and compared to 
equilibrating the system with a finite-temperature bath. 

To our knowledge, this is the first time these features 
have been systematically explored as open-loop control 
problems and solved in a minimal setting by coherent 
local controls and bang-bang modulation of a single local 
noise source that is exactly Markovian. For state transfer, 
our open-loop protocol ensures full state controllability, 
so here it is as powerful as the closed-loop measurement- 
based feedback scheme in Q (see the Supplement |23l . 
IF]). Thus it may serve to simplify many experimental 
implementations . 

We have extended our optimal-control platform DY- 
NAMO [Io| by controls over Markovian noise ampli- 
tudes. Our method was also shown to supersede algo- 
rithmic cooling [U ■ As possible applications, we demon- 
strated the initialisation step of quantum computing (i.e. 
the transfer from the thermal state to the pure zero- 
state |Tl|). the erasure, and the interconversion of ran- 
dom pairs of mixed states, as well as the noise-controlled 
generation of maximally entangled states. 

We anticipate that our approach of switching or even 
modulating the amplitudes of standard Markovian noise 
processes as additional open-loop control parameters (in 
an otherwise coherently controllable system) will pave 
the way to many experimental applications. For instance, 
bit flips may be induced by external random processes 
and amplitude damping may be mimicked by pumping. 
If needed to facilitate experimental implementation, our 
algorithm can be specialised such as to separate dissi- 
pative and unitary evolution. Otherwise, the algorithm 
parallelises coherent and incoherent controls to an extent 
usually going beyond analytical tractability. 
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Appendix A: Proofs of the Main Theorems and Generalisations 

Theorem 1. Let E a be an n-qubit bilinear control system as in Eqn. ([I]) satisfying Eqn. (J4|) for 7 = 0. Suppose the 
n qubit (say) undergoes (non-unital) amplitude-damping relaxation, the noise amplitude of which can be switched 
in time between two values as "f(t) G {0,7*} with 7* > 0. // the free evolution Hamiltonian Hq is diagonal (e.g., 
Ising-ZZ type), and if there are no further sources of decoherence, then the system T, a acts transitively on the set of 
all density operators pOS 1 : 



reacf) Sa (p ) = pos ± for all p € poSj , (Al) 
where the closure is understood as the limit 7*T — > 00. 

Proof. We keep the proofs largely constructive. By unitary controllability, po may be chosen diagonal as 
po =: diag(ro). Since a diagonal po commutes with a diagonal free evolution Hamiltonian Hq, the evolution 
under noise and coupling remains purely diagonal, following 



r(t) 



r =:R a (t)r , (A2) 



where e := e~ 7 *' and R a (t) is by construction a stochastic matrix. With the noise switched off, full unitary control 
includes arbitrary permutations of the diagonal elements. Any of the pairwise relaxative transfers between diagonal 
elements pa and pjj (with i =/= j) lasting a total time of r can be neutralised by permuting pa and pjj after a time 



7; 



1 ln ( P^ +1 ' T + Pjj ) (A3) 



Pjj 



and letting the system evolve under noise again for the remaining time r — Ty. Thus with 2' — 1 such switches all 
but one desired transfer can be neutralised. As p(t) remains diagonal under all permutations, relaxative and coupling 
processes, one can obtain any state of the form 

P (t) = diag(. . . , [ Pll + Pjj ■ (1 - e-"i*% v . . . , [ Pn ■ e-i'%, . . .). (A4) 

Sequences of such transfers between single pairs of eigenvalues pa and pjj and their permutations then generate (for 
7*T — > 00) the entire set of all diagonal density operators A C pos 1 . By unitary controllability one gets all the unitary 
orbits U(A) = pos 1 . Hence the result. □ 

Theorem 2. Let be an n-qubit bilinear control system as in Eqn. |lj satisfying Eqn. ((4]) (7 = 0) now with the n th 
qubit (say) undergoing (unital) bit-flip relaxation with switchable noise amplitude j(t) G {0,7*}. If the free evolution 
Hamiltonian Hq is diagonal (e.g., Ising-ZZ type), and if there are no further sources of decoherence, then in the limit 
"f^T — > 00 the reachable set to £5 explores all density operators majorised by the initial state po, i.e. 



Ccacf) E6 (po) = {p e pos 1 | p -< p } for any p Q e pos 1 . (A5) 
Proof. Again consider the initial state po =: diag(ro). The evolution under the noise remains diagonal following 

j(n-i) _ if(l + e) (1-ef 



r(t) = 



12 ^2l(i_ e ) (i + e ) 



r =: Rb(t)r , (A6) 



where e := e 2 ' and Rb{t) is doubly stochastic. In order to limit the relaxative averaging to the first two eigenvalues, 
first conjugate po with the unitary 

1 -i^ 2 - 1 - 1 



U 12 :=h®^\: 1 j (A7) 

to obtain p' := Ui2PqU\ 2 . Then the relaxation acts as a T-transform [2^| on the first two eigenvalues of p' Q , while 
leaving the remaining ones invariant. 
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Yet the protected subspaces have to be decoupled from the free evolution Hamiltonian Ho assumed diagonal. Any 
such H decomposes as Hq =: i? ,i ® 1 + H 02 ®<t z , where H 01 and H 02 are diagonal. The term with H i commutes 
with p' and can thus be neglected. The other term can be sign-inverted using 7r-pulses in the x-direction on the noisy 
qubit (generated by H\ x ), 

e iK-H lxe -t(r+iH 0:2 ®°-z) e -iTrH 1 x _ e -t(r-iH 0:2 ^ z ) (A8) 

which also leave the bit-flip noise generator invariant. Thus Hq may be fully decoupled in the Trotter limit 

lim ( e -iki T + lfl o.2®<Jz) e -^:(T-tH (K2 ®*^k = e -tr _ ^ A9 ^ 

k— >-oo 

By combining permutations of diagonal elements with selective pairwise averaging by relaxation, any T-transform of 
Po [13 can t>e obtained in the limit — > oo: 

7* 

p{t) = diag ( . . . , \\pu + pjj + (pa - pjj) ■ e~T ] u , . . . , 

7* 

| [pa + Pjj + (Pjj - Pa) ■ e~ % }jj , . . . ) . (Af 0) 

Now recall that a vector y £ M. N majorises a vector x £ R , x ~i y, if and only if there is a doubly stochastic matrix 
D with x = Dy, where D is a product of at most N — 1 such T-transforms (e.g., Thm. B.6 in J28j or Thm. II. 1.10 
in Actually, by the work of Hardy, Littlewood, and Polya [13 this sequence of T-transforms is constructive [HI, 

p32] as will be made use of later. Thus in the limit 7*T — > oo all diagonal vectors r -< ro can be reached and hence 
by unitary controllability all the states p -< po- 

Finally, to see that one cannot go beyond the states majorised by the initial state, observe that controlled unitary 
dynamics combined with bit-flip relaxation is still completely positive, trace-preserving and unital. Thus it takes the 
generalised form of a doubly-stochastic linear map $ in the sense of Thm. 7.1 in [3lj . which for any hermitian matrix A 
ensures $(A) -< A. Hence the (closure of the) reachable set is indeed confined to p -< p . □ 



The conditions for the drift Hamiltonian Hq in theorems above can be relaxed to the following generalisations. 
Any free evolution Hamiltonian Hq may be diagonalised by a unitary transformation: Hq = UHq^U^. The same 
transformation U, when applied to the Lindblad generator V , yields a new Lindblad generator V := UVW . If V 
satisfies Theorem [1] or [2] with any diagonal free evolution Hamiltonian, then V' will satisfy them with Hq. Degenerate 
eigenvalues of Hq yield some freedom in choosing U which, together with arbitrary permutations, can be used to 
make V simpler to implement (e.g. local). 

Moreover, the theorems above are stated under very mild conditions. So 



1. the theorems hold a forteriori if the noise amplitude is not only a bang-bang control jit) £ {0,7*}, but may 
vary in time within the entire interval j(t) £ [0,7*]; the use of this will demonstrated in more complicated 
systems elsewhere; 

2. likewise, if several qubits come with switchable noise of the same type (unital or non-unital), then the (closures 
of the) reachable sets themselves do not alter, yet the control problems can be solved more efficiently; 

3. needless to say, a single switchable non-unital noise process (equivalent to amplitude damping) on top of unital 
ones suffices to make the system act transitively; 

4. for systems with non-unital switchable noise (equivalent to amplitude damping), the (closure of the) reachable 
set under non-Markovian conditions cannot grow, since it already encompasses the entire set of density operators 
(see Sec. [F]) — yet again the control problems may become easier to solve efficiently; 

5. likewise in the unital case, the reachable set does not grow under non- Markovian conditions, since the Markovian 
scenario already explores all interconversions obeying the majorisation condition (see also Sec. IF)) ; 

6. the same arguments hold for a coded logical subspace that is unitarily fully controllable and coupled to a single 
physical qubit undergoing switchable noise. 
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Appendix B: Generalisation of Noise Generators and Their Relation to Coupling to Finite- Temperature 

Baths 

1. Generalised Lindblad Terms 



The noise scenarios of the previous theorems can be generalised by using the Lindblad generator Vg 
with 9 € [0, 1]. Using the short-hand 9 := 1 — 9, the Lindbladian and its exponential turn into 



(1-9) 
9 



T(9) 



I-9 2 




















\ 



and 



(Bl) 



7 



-7.tr(e) 



(cg(9 2 + 9 2 e ) 

e' g cosh(j J9_9) e' g smh(^ 

e' e smh(^*t99) cgcosh^ 
\cg6\l-eg) 



cg9 2 (l-e e )\ 


_ 

cg(9 2 + 9 2 e 9 )J 



(B2) 



with eg := — 1 , eg := e t**/ c " and e' g := e 1 *^ 90 1 I 2 * 1 as further short-hands. Choosing the initial state diagonal, 

9 -{-6 

the action on a diagonal vector of an n-qubit state takes the following form that can be decomposed into a (scaled) 
convex sum of a pure amplitude-damping part and a pure bit-flip part (cp. Eqs. (jA2[) , (|A6[0 



Rg(t) = if ( ™- X) 



C() 



9 2 + 9 2 e e ) 9 2 {l-eg) 
l 2 {l-eg) (9 2 + 9 2 eg) 



1 (I -eg) 
g2 +e 2 \0 eg 



(1 + eg) (I -eg) 
+* 2 Ul-ee) + 



(B3) 



In order to limit the entire dissipative action over some fixed time r to the first two eigenvalues (as in Thin. 1), one 
may switch again as in Eqn. (|A3|) after a time 



eg , I e 
— in 



(e 2 Pv 



l )(Pu + Pjj) 



This is meaningful as long as < Ty(^) < r, which corresponds to the condition 

2 /9 2 < Pii / Pjj < 9 2 /9 2 . 
If 9 ^ 1/2 the noise qubit has a unique fixed point, 

Poo \P) - Cg I Q ^ 2 



(B4) 



(B5) 



(B6) 



Comparing this with the canonical density operator of a qubit with energy level splitting A at inverse tempera- 
ture f3 



k B T> 



1 



=/3A/2 







PP : ~ 2cosh(/3A/2) V e-^' 2 
we can see that the parameter 9 corresponds to the inverse temperature 



(B7) 



P{ff) = | artanh (6(9)) with 6(9) 



9 2 + 9 2 



(B8) 



Thus (for t — > oo) the relaxation by the single Lindblad generator Vg shares the canonical fixed point with equilibrating 
the system via the noisy qubit with a local bath of temperature /3(9). As limiting cases, pure amplitude damping is 
brought about by a bath with zero temperature Tg = (i.e. 9 = 0), while pure bit- flip shares the canonical fixed point 
with the high-temperature limit Tg — > oo (i.e., 9 5), see also Sec. IB 21 for the relation to heat baths. 
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Whereas in a single-qubit system with unitary control and bang-bang switchable noise generator Vg it is straight- 
forward to see that one can (asymptotically) reach all states with purity less or equal to the larger of the purities of 
the initial state po and p oc (6), 



«acf) lqubit:Efl (p ) = {p\p-< Po} U {p'\p < Poo(Q)} , (B9) 

the situation for n > 2 qubits is more involved: relaxation of a diagonal state can only be limited to a single pair of 
eigenvalues, if all the remaining ones can be arranged in pairs each satisfying Eqn. (|B5|) . 

However, Eqn. (|B5j) poses no restriction in an important special case, i.e. the task of cooling: starting from 
the maximally mixed state, optimal control protocols with period-wise relaxation by V$ interspersed with unitary 
permutation of diagonal density operator elements clearly include the partner-pairing approach [24| to algorithmic 
cooling with bias 6(6) defined in Eqn. (|B8[) as long as < 6 < 1/2. Note that this type of algorithmic cooling 
proceeds also just on the diagonal elements of the density operator, but it involves no transfers limited to a single pair 
of eigenvalues. Let ps define the state(s) with highest asymptotic purity achievable by partner-pairing algorithmic 
cooling with bias 6. As the pairing algorithm is just a special case of unitary evolutions plus relaxation brought about 
by Vg, one arrives at 

reacf)(p ) 3 reocf)(p 5 ) for any p , (BIO) 



because any state po can clearly be made diagonal to evolve into a fixed-point state obeying Eqn. (|B5j) , from whence 
the purest state ps can be reached by partner-pairing cooling. 

To see this in more detail, note that a (diagonal) density operator pg of an n-qubit system is in equilibrium with a 
bath of inverse temperature [5(6) coupled to its terminal qubit, if the pairs of consecutive eigenvalues satisfy 

Pl1 = — for all odd i < 2" . (Bll) 

Pi+i,i+i 6 Z 

Hence (for 6 ^ 1/2) such a pg is indeed a fixed point under uncontrolled drift, i.e. relaxation by Vg and evolution 
under a diagonal Hamiltonian Hq thus extending Eqn. IB6I to n qubits. Now, if (say) the first pair of eigenvalues 
is inverted by a selective 7r pulse (which can readily be realized by unitary controls with relaxation switched off), a 
subsequent evolution under the drift term only affects the first pair of eigenvalues as 



1 



6 2 + 6 2 



6 2 + 6 2 e e ) 6 2 (l-e e ) 
I 2 (I -eg) (6 2 + 6 2 eg) 



Cg6 2 \ 1 (9 2 + £g(6 2 - 6 2 ) 

cgd 2 ~ 9 2 + 6 2 W-eg(6 2 - 6 2 ) 



eg (l-£fl) 
(I -eg) eg 



(B12) 

l2\ 



cgd 2 



In other words, the evolution then acts as a T-transform on the first eigenvalue pair. Since the switching condition 
Eqn. IB5I is fulfilled at any time, all T-transformations with eg £ [0, 1] on the first pair of eigenvalues can be obtained 
and preserved during transformations on subsequent eigenvalue pairs. 

Hence from any diagonal fixed-point state pg (including ps as a special case), those other diagonal states (and 
their unitary orbits) can be reached that arise by pairwise T-transforms as long as the remaining eigenvalues can be 
arranged such as to fulfill the stopping condition Eqn. IB5I Suffice this to elucidate why for n > 2 a fully detailed 
determination of the asymptotic reachable set in the case of unitary control plus a single switchable Vg on one qubit 
appears involved and will thus be treated elsewhere. 

2. Bosonic or Fermionic Heat Baths 

Following the lines of [32l , Ch. 3.4.2], a qubit with the Hamiltonian 

H = -^ y (B13) 

coupled to a bosonic (or fermionic) heat bath with inverse temperature (3 via a coupling of the form a x ® (a + a') is 
described in the Born-Markov approximation by a Lindblad equation (of the form of Eqs. JTJ) , with two dissipator 
terms, r a - and T a +, 

P=-(iH + l(l±n(u ))T a - + 7 n( Wo )r CT+ )p, (B14) 
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where n(u>) := l/(e l3hw =p 1) is the Planck (or Fermi) distribution and 7 is the relaxation rate constant. In the 
vec-superoperator representation this yields a Liouvillian which takes the form 



/0 1 \ 



Ft oc 



-i 



-i 

\o -ly 



n(u) ) 



H 








± M 





i±i 

2 













i±i 

2 





V 1 








TV 



In the zero-temperature limit, n(u>o) — > and only the T a - term remains. In the limit T 
for bosons and n(wo) ~> \ for fermions, and one thus obtains in both cases 



lim T T cx r {a + -1 



-10 1 
0-100 
0-10 
10 0-1, 



(B15) 

00, we have ti(ujq) — > 00 
(B16) 



This is equivalent to dissipation under the two Lindblad operators {cr m ,o- y } since ^{a x , a } — ^{a+.a-}- ln contrast, 
o~ x as the only Lindblad operator (generating bit-flip noise) gives 



r 




(B17) 



in agreement with Eqn. (|B1|) . Note, however, that the propagators e~ tT " generated by Ti a +^-\ vs rr^i only act on 
diagonal density operators (and those with purely imaginary coherence terms p\2 = P21) m an indistinguishable way. 
However, the relaxation of the real parts of the coherence terms pi2 = p\\ differs: Y am leaves them invariant, see also 



Eqn. flB2J 
while Too 
T(6). 



while r 



T -j does not. 



In other words, r CT;r does have nontrivial invariant subspaces used in Eqn. (|A7 



does not. Also for non-zero temperatures, there is no direct equivalence between relaxation under Ft and 



Appendix C: grape Extended by Noise Controls 

In state transfer problems the fidelity error function used in [lfj | is valid if the purity remains constant or the target 
state is pure. In contrast to closed systems, in open ones these conditions need not hold. Thus here we use a full 
Frobenius-norm distance-based error function instead: Sp := ||Am : o — -^taxgetllp i where Xk-.o = Xk ■ ■ ■ Aivec(po) is 
the vectorised state after time slice k, Xk = e~ AtLk is the propagator for time slice k in the Liouville space, and 
Lfc := iH u (tk) + r(£fc). The gradient of the error is obtained as 

- - 2Retr f (X M:0 - X^)*™™) , (CI) 



auj(tk) V auj(tk 

where 

OXm-.q dX k 

— — = X M ■ ■ ■ X fe+ i-— — -X fe _i • • • Xivec (po . C2 

OUj{t k ) OUjitk) 

Exact expressions for the derivatives of Xk [T3| require L k to be normal, which does not hold in the general case of 
open systems of interest here. Instead we may use, e.g., the finite- difference method to compute the partial derivatives 

as 

dX k _ exp(-Ai(L fe + isHj)) - exp(-AtL fc ) 



lim ~ yv v K ' 311 ^ K > , (C3) 

dX k exp(-A<(L fe + s r Vj ))-exp(-AtL fc ) 

- lim . (C4) 



9lj(tk) s^o 

The optimal value of s is obtained as a tradeoff between the precision of the gradient and numerical accuracy, which 
starts to deteriorate when s becomes very small. 
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(a) (b) (c) 




Figure 4. Same as in Fig. [T] but with additional non-switchable background dephasing noise on all the three qubits. (a) Quality 
versus dephasing rate 70 (with 70/ J £ {0.01, 0.02, 0.05, 0.1, 0.2}) for sequences of duration T — 8/J. The dots (red circles for 
averages) are individual numerical optimal-control runs with random initial sequences, (b) Evolution of the eigenvalues under 
the best sequence for the strongest background noise (70 = 0.2J) leading to the zero-state with a considerably low residual 
error of Sf — 0.077. This sequence (c) shows five relaxative periods with maximal noise amplitude on qubit one (701) for 
transforming eigenvalues, while the unitary actions again mainly take place in the intervals between them. 



Appendix D: Further Numerical Results 

Example lb. Interestingly, the initialisation task of Example 1 can still be accomplished to a good approximation 
when unavoidable constant dephasing noise on all the three qubits is added. This is shown in Fig. 3] for a range of 
dephasing rate constants reaching from 1% to 20% of the coupling constant. Though the dephasing does not affect 
the evolution of diagonal states, it interferes with the i-swaps needed to permute the eigenvalues. For 7/3 = 0.2 J, 
numerical optimal control suggests the sequence Fig. [4] (c) with five dissipative steps and increasing time intervals for 
the z-swaps. 

In Example 2, we considered erasing the pure initial state P|oo...o) to the thermal state pth by controlled bit-flip 
noise of Eqn. ([3]) to illustrate the scenario of Thin. [5] For n qubits, one may use a similar n-step protocol as in 
Example 1, this time approximately erasing each qubit to a state proportional to 1. Again one finds that the residual 
error 5f is minimal for equal T q to give 

^) = ^((i + e T-i), (Dl) 

where e := e -7.T„/(2n)_ Thig yieldg 

T b = Q j - £ In ((2 n S 2 Fb + l) 1 /" - 1) . (D2) 

Once again Fig. OJa) shows that numerical optimal control finds much faster solutions than this simplistic protocol. 
The noise amplitude tends to be maximised throughout the sequence with the unitaries fully parallelised, as shown in 
the example sequence (c), and reflected in the eigenvalue flow (b). This works so well because p t h is the unique state 
majorised by every other state, and thus all admissible eigenvalue transfers lead towards the goal. 

The advantage of optimal-control based erasure becomes evident when comparing it to free evolution: Pure bit- 
flip noise on one qubit (without coherent controls) would just average pairs of eigenvalues once if the free-evolution 
Hamiltonian is mere Ising-ZZ coupling, which commutes with the initial state. Hence free evolution does not come 
closer to the thermal state than 8f — 0.61 and only by allowing for unitary control, erasure becomes feasible for the 
Ising chain. 



Appendix E: Analytical Scheme Following Hardy, Littlewood, and Polya 

The work of Hardy, Littlewood, and Polya [30j | (HLP) provides a constructive scheme ensuring the majorisation 
condition /Otargct ~< p(t) -< Po f° r all < t < T to be fullfilled for all intermediate steps. Let the initial and the 
target state be given as diagonal vectors with the eigenvalues of the respective density operator in descending order, 
so po =: diag(yi, j/2 7 • ■ • 7 2/tv) and ptarget =: diag(xi, 22, ■ • ■ , xn)- Following [28|, p32f], fix j to be the largest index 
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(a) (b) (c) 




sequence duration T [1/J] time[1/J] time[1/J] 



Figure 5. Transfer from the zero-state P|ooo) to the thermal state pth = |l in a 3-qubit Ysmg-ZZ chain with controlled bit-flip 
on qubit one and local a;,j/-pulse controls on all qubits as in Example 2. (a) Quality versus total duration T, with the dashed 
line as the upper bound from Eqn. (|D2[) . Dots (red circles for averages) denote individual numerical optimal-control runs 
with random initial sequences, (b) Evolution of the eigenvalues under the controls of the best of the T = 3/ J solutions. The 
corresponding control sequence (c) shows that the noise is always maximised, and the unitary actions generated by (u xv ,u yv ) 
are fully parallelised with it. 



such that Xj < yj and let k > j be the smallest index with x^ > y^. Define r5 := min{(i/j — Xj), (x^ — yk)} and 
A := 1 — 8/{y,j — yk). This suffices to construct 

y' ■■= Mj + (i - X)Qjk v (El) 

satisfying x -< y' -< y. Here the pair-permutation Qjk interchanges the coordinates y^ and yj in y. So y' is a T- 
transform of y, and Ref. [28[ shows that by N — 1 successive steps of T-transforming and sorting, y is converted into 
x. Now the T-transforms A1+ (1 — \)Qjk can actually be brought about by switching on the bit-flip noise according 
to Eqn. IA10I for a time interval of duration 

r jfc :=-|-ln|l-2A|. (E2) 

With these stipulations one obtains an iterative analytical scheme for transferring any p by unitary control and 
switchable bit-flip noise on a terminal qubit into any ptarget satisfying the reachability condition ptarget -< Po- 



Scheme for Transferring Any ri-Qubit Initial State po into Any Target State ptargct -< po by Unitary 
Control and Switchable Bit-Flip Noise on Terminal Qubit: 

(0) switch off noise to 7 = 0, diagonalise target U x ptargetC/J =: diag(a;) to obtain diagonal vector in descending order 
x — (xi,X2, ■ ■ ■ ,xn); keep U x ; 

(1) apply unitary evolution to diagonalise po and set po =: diag(y); 

(2) apply unitary evolution to sort diag(y) in descending order y = (2/1,3/2, ■ ■ ■ , yjv); 

(3) determine index pair (j, k) by the HLP scheme; 

(4) apply unitary evolution to permute entries (j/i, yj) and (j/2, J/fe) of y, so diag(y) = diag(j/j, yk, ■■■)', 

(5) apply unitary evolution U12 of Eqn. (|A7[) to turn p y — diag(y) into protected state; 

(6) switch on bit-flip noise on terminal qubit j(t) = 7, for duration Tjk of Eqn. ()E2|I (while decoupling as in 
Eqn. (|A9|l 'l to obtain p y r, 

(7) to undo step (5), apply inverse unitary evolution U{ 2 to re-diagonalise p y i and obtain next iteration of diagonal 
vector y = y' and p y = diag(y); 

(8) go to (2) and terminate after N - 1 loops (N := 2™); 

(9) apply inverse unitary evolution Jj\ from step (0) to take final p y to U\p y U x ~ ptargct- 

Note that the general HLP scheme need not always be time-optimal: E.g., a model calculation shows that just the 
dissipative intervals for transferring diag(l, 2, 3, ... , 8)/36 into lg/8 under a bit-flip relaxation-rate constant 7* = 5J 
and achieving the target with Sp = 9.95 • 10 -5 sum up to T rc i ax = 12/ J in the HLP-scheme, while a greedy alternative 
can make it within T r ' olax = 6.4/J and a residual error of 6f = 6.04 ■ 10 -5 . 
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Appendix F: Outlook on the Relation to Extended Notions of Controllability in Open Quantum Systems 

The current results also pave the way to an outlook on controllability aspects of open quantum s yst ems on a more 
general scale, since they are much more intricate than in the case of closed systems [3, 13.171 Il5l. [iH l2ll I22I, |33T-[35| . 

Here we have taken profit from the fact that like in closed systems (where pure-state controllability is strictly weaker 
than full unitary controllability [36l |37|). in open quantum systems Markovian state transfer appears less demanding 
than the operator lift to the most general scenario of arbitrary quantum map generation (including non-Markovian 
ones) first connected to closed-loop feedback control in Therefore in view of experimental implementation, the 
question arises how far one can get with open-loop control including noise modulation and whether the border to 
closed-loop feedback control is drawn by Markovianity. 

Due to their divisibility properties [H, [3§| that allow for an exponential construction (of the connected component) 
as Lie semigroup [l6| . Markovian quantum maps are a well-defined special case of the more general completely positive 
trace-preserving (CPTP) semigroup of Kraus maps, which clearly comprise non-Markovian ones, too. While some 
controllability properties of general Kraus-map generation have been studied in [H, [||, a full account of controllability 
notions in open systems should also encompass state-transfer to give the following major scenarios: 

1. Markovian state-transfer controllability (MSC), 

2. Markovian map controllability (MMC), 

3. general (Kraus-map mediated) state controllability (KSC) (including the infinite-time limit of 'dynamic state 
controllability' (DSC) @,0|), 

4. general Kraus-map controllability (KMC) 

Writing 'C' and 'Q in some abuse of language for 'weaker than' and 'strictly weaker than', one obviously has at least 
MSC C KSC and MMC C KMC, while DSC S KMC was already noted in the context of control directly over the Kraus 
operators Q. In pursuing control over environmental degrees of freedom, Pechen 0, |40| also proposed a scheme, 
where both coherent plus incoherent light (the latter with an extensive series of spectral densities depending on ratios 
over the difference of eigenvalues of the density operators to be transferred) were shown to suffice for interconverting 
arbitrary states with non-degenerate eigenvalues in their density-operator representations. 

Yet the situation outlined above is more subtle, since unital and non-unital cases may differ. In this work, we have 
embarked on unital and non-unital Markovian state controllability, MSC (4l| : 

Somewhat surprisingly, in the non-unital case (equivalent to amplitude damping), the utterly mild conditions of 
unitary controllability plus bang-bang switchable noise amplitude on one single internal qubit (no ancilla) suffice for 
acting transitively on the set of all density operators (Theorem [1}. Hence these features fulfill the maximal condition 
KSC already. In other words, for cases of non-unital noise equivalent to amplitude damping (henceforth indexed by 
'nu'), KSC nu implies KSC. Moreover, under the reasonable assumption that the mild conditions in Theorem Q] are in 
fact the weakest for controlling Markovian state transfer MSC nu in our context, Theorem Q] shows that MSC nu implies 
KSC via KSC nu . So in the (extreme) non-unital cases, there is no difference between Markovian and non-Markovian 
state controllability. — On the other hand in order to compare non-unital with unital processes, taking Theorems [T] 
and [2] together proves MSC U S MSC nu , since the former is restricted by the majorisation condition of Theorem^ 

Similarly, in the unital case (equivalent to bit- flip), the mild conditions of unitary controllability plus bang-bang 
switchable noise amplitude on one single internal qubit suffice for achieving all state transfers obeying majorisation 
(Theorem [5]) • Hence again they fulfill the maximal condition KSC U at the same time. This is because state transfer 
under every unital CPTP Kraus map (be it Markovian or non-Markovian) has to meet the majorisation condition; so 
we get KSC U . On the other hand, the majorisation condition itself imposes the restriction KSC U £ KSC. Again, under 
the reasonable assumption that the mild conditions in Theorem [5] are in fact the weakest for controlling Markovian 
state transfer MSC U in our context, Theorem [5] shows that MSC U implies KSC U . Thus also in the unital case, there is 
no difference between Markovian and non-Markovian state controllability. 

The results on these two cases, i.e. non-unital and unital (in the light of Appendix IBl seen as the limits 6 = and 
9 = h, respectively), can therefore be summarized as follows: 

Corollary 1. In the two scenarios of Theorem\l\ (non-unital) and\^ (unital), Markovian state controllability already 
implies Kraus-map mediated state controllability and one finds 

MSC nu =^> KSC nu =^> KSC 

Utf (Fl) 
MSC U =>> KSQ, 
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However, whether MSQ =>■ KSQ also holds in the generalisation of Appendix iBl where 9 can range over the entire 
interval 9 € [0, ^} (with 9 = giving the limiting cases MSC nu , KSC nu and 9 = \ yielding MSC U , KSC U ), currently 
remains an open question. 

This has an important consequence for experimental implementation of state transfer in open quantum systems: On 
a general scale in ?i-qubit systems, unitary control plus measurement-based closed-loop feedback from one resettable 
ancilla (as, e.g., in Ref. [H following Q) can be replaced by unitary control plus open- loop bang- bang switchable 
non-unital noise (equivalent to amplitude damping) on a single internal qubit. This is because both scenarios are 
sufficient to ensure Markovian and non-Markovian state controllability KSC. Example 5 in the main part illustrates 
this general simplifying feature. 

Yet some questions with regard to the operator lift to map synthesis remain open: Assessing a demarcation between 
MMC and KMC (and their unital versus non-unital variants) seems to require different proof techniques than used here. 
In a follow-up study we will therefore further develop our lines of assessing the differential geometry of Lie semigroups 
in terms and their Lie wed ges 
known to be more intricate 139 
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22j to this end, since judging upon Markovianity on the level of Kraus maps is 
42j. More precisely, time- dependent Markovian channels come with a general form of 
a Lie wedge in contrast to time-independent Markovian channels, whose generators form the special structure of a Lie 
semialgebra (i.e. a Lie wedge closed under Baker-Campbell-Hausdorff multiplication). In [ll|, we have therefore drawn 
a detailed connection between these differential properties of Lie se mig roups and the different notions of divisibility 
studied as a defining property of Markovianity in the seminal work [38j . 

Again, these distinctions will decide on simplest experimental implementations in the sense that measurement- 
based closed-loop feedback control may be required for non-Markovian maps in KMC, while open-loop noise-extended 
control may suffice for Markovian maps in MMC. More precisely, closed-loop feedback control was already shown to 
be sufficient for KMC in Q (which was the aim that work set out for), yet it remains to be seen whether it is also 
necessary, and in particular, if it is necessary for the (supposedly) weaker notion MMC. If it turns out not to be 
necessary, then measurement-based closed-loop feedback control on a system extended by one resettable ancilla 
would be not be stronger than our open-loop scenario of full unitary control extended by (non-unital) noise modulation 
not only in the case of state transfer, but also in quantum-map synthesis with direct bearing on the simplification of 
quantum simulation experiments Q. 



